PRODUCT VACUA WITH BOUNDARY STATES AND THE 
CLASSIFICATION OF GAPPED PHASES 



SVEN BACHMANN AND BRUNO NACHTERGAELE 

Abstract. We address the question of the classification of gapped ground states in one dimension 
that cannot be distinguished by a local order parameter. We introduce a family of quantum spin 
systems on the one-dimensional chain that have a unique gapped ground state in the thermody- 
namic limit that is a simple product state but which on the left and right half-infinite chains, have 
additional zero energy edge states. The models, which we call Product Vacua with Boundary States 
(PVBS), form phases that depend only on two integers corresponding to the number of edge states 
at each boundary. They can serve as representatives of equivalence classes of such gapped ground 
states phases and we show how the AKLT model and its SO(2J + l)-invariant generalizations fit 
into this classification. 



1. Introduction and statement of the results 

In this paper we apply the notion of automorphic equivalence introduced in p] to a class of 
quantum spin chains with nearest neighbor interactions. We view this as a step toward a qualitative 
theory of gapped ground state phases of quantum many-body systems, with the ultimate goal of 
understanding these phases and the transitions between them as well as we do phase transitions 
at finite temperature. One significant difference is the occurrence of so-called topological phases, 
i.e. , the situation where different ground states coexists that cannot be distinguished by a local 
observable (order parameter). This situation is in stark contrast to the classical case of spontaneous 
symmetry breaking for which an order parameter distinguishing the different pure phases is readily 
identified. 

The topological degeneracy of the ground states in models where it occurs, is usually best under- 
stood by considering the models defined on a set of lattices, rather than just one lattice. Members 
of this set of lattices can have boundaries or non-trivial topology. The discovery about a decade 
ago [2] of two-dimensional systems with topological order sparked great interest, specifically in view 
of possible applications in quantum information processing. Since then, however, this interest has 
broadened to aim more generally at understanding the structure of gapped ground state phases 
and the quantum phase transitions between them. 

The very criterion that determines whether two models should be considered to be in the same 
ground state phase has only recently received careful attention, both from a physical point of 
view [3] and from a more mathematical one [U [5]. For systems that have a spectral gap in the 
thermodynamic limit, we have proposed in [l] the notion of local automorphic equivalence to define 
a phase and shown how it relates to other criteria. Simply put, two models are in the same gapped 
ground state phase if there exists a quasi-local automorphism of the algebra of observables mapping 
the ground state space of one of the system to that of the other. Here, quasi-local is meant in the 
strong sense that the automorphism satisfies a Lieb-Robinson bound with error estimates that 
vanish faster than any power law. 

The notion of automorphic equivalence becomes distinct from unitary equivalence, and more 
useful, when one applies it to infinitely extended systems. To study topological order one then 
also considers Hamiltonians not only on one lattice, but on a family of lattices with boundaries 
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or non-trivial topologies and the automorphisms map sets of ground states onto each other. This 
implies in particular that for any given lattice the ground state degeneracies are equal for models 
that are in the same phase. 

This work is concerned with one-dimensional spin chains. The different types one-dimensional 
lattices we consider are (i) finite chains, (ii) left and right half-infinite chains, (iii) and the full 
bi-infinite chain. Despite the absence of interesting topology, this one dimensional situation allows 
us to emphasize the role played by boundaries and edge ground states. This yields a physically 
relevant refinement of the consensus that all gapped one-dimensional models with a unique ground 
state in the thermodynamic limit are in the same phase [3], usually referred to as the Haldane 
phase. We will show that models that share the same and simplest state on the doubly infinite 
chain can differ by their boundary behavior on the two possible half- infinite chains. These examples 
support the view that such systems can be classified by their boundary behavior. 

To make this concrete, let us consider Hamiltonians H\ defined on finite subsets of Z by 

XcA 

where, for each finite subset of Z, $(X) is a self-adjoint element of the algebra of observables 
belonging to X, which will be denoted by Ax- The dynamics and ground states of the model on 
the infinite lattice Z, and the half-infinite lattices (—00, 0] and [1, +00), are defined on the algebras 
of quasi-local observables Az, «4.(-oo,01) an d *4[i,+oo)j respectively. 

Let {T x } x be the group of translations on Z and {r^j^ be the representation of the translations 
as automorphisms of the quasi- local algebra of observables Az- We assume that $ has finite range 
and is translation invariant, <&(T X (X)) = t x ($(X)), for all X C Z. For the infinite systems on the 
lattices V = Z, (—00, 0], [1, +00), a ground state is a state of the quasi-local algebra „4r such that 

uj{X*[H,X}) > 

for all local observables X G Ar- Let S r be the set of ground states on T. All states obtained 
as thermodynamic limits of finite- volume ground states are ground states of the infinite system in 
this sense. In order not to unnecessarily complicate things, we will say that the system is gapped 
if there exists 7 > such that the difference between the smallest and the next-smallest eigenvalue 
of is bounded below by 7, for all finite intervals A. We will denote by S^, the set of states 
of a quantum system with a d-dimensional Hilbert space. E.g., £2 is the standard Bloch sphere 
describe all possible states of a qubit, and £\ contains a single pure state. For two convex sets of 
states A and B, A = B will mean isomorphism of convex sets, and for two sets of states Sq and S\ 
of a given algebra of observables A, Sq ~ S\ will mean 

(1.1) 5i = {wo a |wG Sq}, 

with a a quasi-local automorphism of A, satisfying 

||[a(A),S]|| < \\A\\ (d(supp(A),supp(B)), 

for a function F that decays faster than any power law, d is the standard distance between subsets 
if the lattices and supp(C) denotes the smallest set X, such that C G Ax- 

Conjecture 1.1. Let <3?o and 3>i be two translation invariant finite range interactions. Assume 
that for i = 0, 1, there exist integers k and r{ such that 

Then, the following are equivalent: 

i. r = n,Zo = h; 

ii. There exists a continuously differentiate family of translation invariant finite range inter- 
actions <3?(s), with s £ [0, 1], such that $(0) = $0 an d ^(1) = $1, and there is a uniform 
lower bound 7 > for the gap of all models defined by <3?(s), s G [0, 1]. 
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Clearly, (iii) implies (i), since the automorphism a in (1.1) induces a homeomorphism between 
two sets of states of the form regarded as topological manifolds. That (ii) implies (iii) is proved 
constructively in [lj. The remaining implication is conjectural. The condition that the sets of 
ground states 5 r are isomorphic to for some positive integer d, is equivalent to assuming that 
they are the state spaces of the full matrix algebra of complex d x d matrices. It seems plausible 
that the assumption that 5 Z is a single point, which means a unique ground state in the bulk, 
implies this condition. In this paper, we limit ourselves to this situation, but it is straightforward 
to generalize the conjecture to cases where the sets 5 r are isomorphic to the state spaces of a 
direct sum of matrix algebras. E.g., if the bulk ground states is two-fold degenerate, i.e. , two 
pure ground states in the bulk, S z is isomorphic to the segment [0,1], which is the state space 
of C C. We then expect that 5* -00 ' ' and 5t 1 ' +0 °-', if finite-dimensional, will be isomorphic to 
the state spaces of the direct sum of two square matrix algebras. This is the situation needed to 
describe the standard (non-topological) ground state phase transitions. 

The duality between 'topological order' in the bulk and the existence of non trivial edge states is 
expected to extend to two-dimensional systems. See for example [HE] for a detailed study of a class 
of examples. In fact, the existence of gapless edge modes in a system that is gapped in the bulk was 
originally proposed as a way to measure topological order in fractional quantum Hall systems [8]. 
Further, the equality of bulk and edge conductances in the two-dimensional quantum Hall effect [9], 
or the bulk-edge correspondence for topological insulators [10] can be seen as instances of the same 
principle. 

We start this paper by presenting a detailed study of the so-called PVBS models we introduced 
in [TT]. These are a continuous family of translation invariant Hamiltonians with nearest neighbor 
interaction and finitely correlated ground states. The spectral gap does not close in the thermody- 
namic limit except when specified parameters take on a critical value (which we have normalized 
to be 1). On the infinite chain, they share a unique product ground state. However, they have 
degenerate ground states on half-infinite chains which have a simple interpretation. There are n 
distinguishable particles that can be added to the vacuum without raising the energy if no other 
particle of the same type is already present and if they can bind to a boundary. Among the n 
types of particles, ni can bind to a left edge while ur = n — ni bind to a right edge, yielding a 
degeneracy of 2 n on finite chains, and 2 nL , resp 2 nR , on the two possible half-infinite chains. We 
call these models Product Vacua with Boundary States. We illustrate the classification of phases 
by proving that two PVBS models belong to the same phase if and only if they share the same 
number of left- and right- binding particles. The simplicity of the models make them good candi- 
dates as representatives of the gapped phases and it follows from our results that Conjecture |1.1| 
holds within this restricted class of models. 

Next, we show that the well-known spin-1 antiferromagnetic chain known as the AKLT model |12j 
fits in the PVBS classification, namely that it belongs to the PVBS phase with one of each boundary 
particle. These proofs rely on two ingredients: First that a local automorphism can be constructed 
from a smooth path of gapped Hamiltonians interpolating between the two models in question. 
And second on a criterion for models with matrix product ground states (MPS) to have a spectral 
gap [13] . Despite the fact that much can be said by simply considering the algebraic deformation 
of the matrices generating the ground state spaces as in [5J, we shall emphasize that this may not 
be sufficient. First of all, the same spaces can be ground state spaces of both gapped and gapless 
frustration- free Hamiltonians, see |14] . Secondly, the pointwise construction of Hamiltonians from 
a smooth family of matrix product maps is not a continuous operation. More fundamentally, 
smoothly connecting MPS described by matrices of different dimensions requires the use of non 
minimal representations in which the spectral gap is in general not under good control. The 
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central objects of the smooth deformation must be the ground state spaces and the Hamiltonians 
themselves. 

Already from the few examples that we study here in detail, it becomes clear that the classification 
of gapped phases that we obtain is finer than those previously pursued in the literature (see, e.g., 
|15l 116].) By this we do not mean that there are more distinct gapped phases in the bulk than 
previously claimed, but our finer analysis provides more information about the ground state phase 
diagrams. To illustrate this, consider two models defined by finite-range interactions that have a 
unique ground state in the bulk and a non-vanishing spectral gap above them. If the half-infinite 
chains with these interactions have inequivalent ground state spaces in the sense of above, any 
smooth curve of finite-range interactions connecting the two given models will contain a point 
where the spectral gap closes. This does not contradict the well-established fact that each model 
can be smoothly connected to a model with a simple product ground state. The interactions of the 
models with the unique product ground state will however not be the same and will, in particular, 
have different edge state behavior. 

One possible family of models generalizing the AKLT model carries the fundamental represen- 
tation of SO{d) and its nearest-neighbor interaction is «SO(d)-invariant [17 . In the last section 
of this article, we propose a smooth path of Hamiltonians between these SO{d) models for odd 
d = 2 J + 1 and the PVBS models with J particles bound to each edge. Here again, the starting 
point is an algebraic deformation, namely of the Clifford algebra underlying the SO(d) invariance 
of the models. 

Other generalizations in one dimension include PVBS models allowing more than one particle 
of each type, generating arbitrary degeneracies at each edge with a unique bulk state, or multiple 
bulk states, namely systems with a broken discrete symmetry. Local continuous symmetries can 
also be included by identifying a subset of particles. Of course, adding symmetries enriches the 
classification of gapped phases and the possibilities of transition between them |15l [16] . 

The class of PVBS models is rich enough to cover a large range of behaviors of one-dimensional 
quantum spin systems. Moreover, their simplicity allows for a complete understanding of their 
ground state properties, and they appear amenable for a general study of the transitions between 
the ground state phases they represent: A phase transition arises in a PVBS model whenever one 
of the edge particles unbinds from one boundary and binds to the other. At the critical point the 
particle becomes a massless excitation of the vacuum and the gap closes. 



2. Product vacua with boundary states 

We consider a quantum spin system on the chain Z. Any point i£Z carries a quantum system 
with Hilbert space C d , for some d G N, equipped with an orthogonal basis {eo,...,e n }, where 
n = d — 1. We shall denote the finite subset [a, b] n Z simply by [a, b], for a < b G Z. The Hilbert 
space %n of a finite chain of length N = b — a + 1 is given by 

H N = ® b x=a C d - 

We first define a nearest-neighbor interaction h G H2 as follows. For < i < j < n, let 

(2.1) ft! = a® ej - (''''• A, 'A, ej <g> Bi , 

(2.2) 4> u = ei®ei, 

where Oij G M, Q%j = —Oji, and / A; £ C. By redefining the phases 6ij, we can assume Aj > 
without loss of generality. The interaction h is defined as the orthogonal projection on the space 
spanned by these vectors. 
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Given a fixed set of parameters A = {Aj:0<i<n} and phases = {6ij : < i < j < n}, the 
PVBS Hamiltonian on [a, b] is given by 

6-1 

H[ at b] = ^2 h x>x+ i 

x=a 

where h x ^ x+ \ is a copy of h acting on the neighboring sites x and x + 1. We consider four different 
realizations of the PVBS models, namely when defined on a finite chain [a, b], on either of the two 
half infinite chains with one edge, and on the infinite chain Z. Of course, the various infinite limits 
involved here must be taken carefully. 

Let 7(A, 0, N) be the spectral gap above the ground state energy of the PVBS model on a chain 
of length N, with parameters A and 0. In the sequel, we shall set Ao = 1. The special role played 
by Ao will be discussed later. 

Proposition 2.1 (Gap bounds). The spectral gap for the Hamiltonian H[ a b j is bounded above by 

(2.3) 7 (A, 9, AO < min { (l - (l + : 1 < i < n 

for any N = b — a + 1, where = (1 + Aj) / 1 1 — Aj| . Moreover, if Aj ^ 1 for all i ^ 0, then 

7 := inf 7 (A,0,AO > 0. 

By standard results [18], the second statement implies a spectral gap for the Hamiltonian in the 
GNS representation of the quasi-local algebra induced by the ground states in the thermodynamic 
limit. The notion of a ground state phase of quantum spin systems was discussed in [TJ, where 
automorphic equivalence between ground state spaces was introduced. Given a differentiable curve 
of local Hamiltonians H(s), for s £ [0, 1], such that there is a lower bound for the spectral gap 
which is uniform in the size of the system and in the parameter s, then there exists a quasi-local 
spectral flow a s y of the algebra of observables mapping the ground state space of H(s') to that of 
H(s). For finite chains it is given by unitary conjugation, and can be extended to infinite systems 
using its good locality properties. 

In the present case, let A^ := {i : Aj < 1}, A^ := {i : Xi > 1}, and 

nL ■= |A L | , n R := \A R \ . 
In the following, we say that two quantum spin systems are equivalent if they satisfy any of the 



condition of Conjecture (1.1). 



Theorem 2.2 (PVBS classes). Two PVBS models with respective are equivalent if and only 
if 



n\ = n\ and n\ = n\. 



Before we enter into the proof of this theorem, let us briefly discuss its importance. The family of 



PVBS models is continuously parametrized by A and 0. Theorem 2.2 distinguishes among them a 
countable set of equivalence classes characterized by a simple and explicit criterion. In fact, all these 
translation invariant models have the same, unique, gapped ground state in the thermodynamic 
limit. However, they do not belong to the same ground state phase. As we will discuss in more 
details later, their properties differ by their boundary behavior which is precisely classified by the 
two integers ul and ur. 

Automorphic equivalence of the models - or equivalently their belonging to the same ground 
state phase - is obtained by constructing a smooth family of gapped Hamiltonians interpolating 
between the two models under consideration, see p]. Reciprocally, the theorem implies that along 
any curve relating two models belonging to two phases indexed by different nL,nji, the gap must 
close and there is a quantum phase transition. 
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The rest of the section is organized as follows. First, we construct the ground states of the 
PVBS models on finite chains. These models are frustration free and finitely correlated, so that 
it will be convenient to use their representation as matrix product states (MPS). The first part of 
Proposition 2.1 follows by a variational argument. Its second statement uses the martingale method 
of [13] and a crucial geometric estimate associated with ground state projections. We construct 



gapped paths of Hamiltonians within the set of PVBS models and Theorem 2.2 follows from the 
gap estimates of Proposition |2.1| 



2.1. Matrix product states (MPS). We briefly review the basic construction and properties 
of purely generated matrix product states, also known as finitely correlated states, see |18j . Let 
e G Aik be a positive element, p £ Ai* k be a state on and let E : Aid ® Aik -M-k De a 
completely positive map such that 

E(l®e) = e, p(E(i ® B)) = p(B) , 

for all B € Aik- The matrix product state uj on the chain algebra A generated by (E, e,p) is the 
extension of 

u{A a ®---®A h ):= p{e)' 1 p (E^ o • • ■ o E Ab (e)) , 

where Ai £ Aid and we have denoted E A (B) = K(A ® B). If E is a pure map, then there exists 
V : C k -> C d ® C k such that E(A ® B) = V* (A ® B)V . 

Now, the breaking or not of translation symmetry can be expressed as a spectral property of the 
transfer operator E := Ej. The state to cannot be decomposed into periodic components iff 1 is an 
isolated, non degenerate eigenvalue of E and all other eigenvalues have a strictly smaller modulus 
- the peripheral spectrum is trivial. This implies the following exponential convergence 

lim EW(B) = &°°\B) := p(B)e 

and since 

u{A a 1 ■ ■ ■ ® 1 ® At) = p{e)~ l p (E Aa o E^- 2 ) o E Afc (e)) , 

the exponential decay of correlations in the ground state. We shall say that the representation of 
uj is minimal if both p and e are invertible. 

The state uj defined on the infinite chain is the thermodynamic limit of 'valence bond vectors' 
defined as follows. The map V can be described by the set {vi}f =1 of k x k matrices such that 

d 

Vx = ^2 ei ® ViX ■ 

1=1 



The range Q N of the map T N :M k -> (C d ) m defined by 



d 



(2.4) T N {B):= Y, Tr ( Bv i 



,»jv=1 



is the set of matrix product states on a finite chain of length N. If {Xp}p=i is an orthonormal basis 
of Aik, a simple calculation yields 

n n 

(T N (B L ),AT N (B R ))= Ti{B L v iN --- Vil )Tr(B R v jN -- ■ Vj 1 ) {e^ ■ ■ ■ ei N , Aej 1 ■ ■ ■ ej N ) 

ii...ijv=Oj'i.-ijv=0 
k 

(2-5) = Y {xpM a N) (Bl \ XP ) (X q \ B R ) Xq ) ■ 

p,q=l 
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In a minimal representation, the spaces Gn always satisfy a crucial intersection property. The 
following lemma presents a sufficient condition for the property to hold which does not require 
faithfulness of the state p. 

We consider the case where the matrices of the MPS are generators of a quadratic algebra. Let 
U := span{<7j : 1 < i < d} and (B^LqU® 01 the tensor algebra generated by IA. Given a multiindex 
/ = . . . ,i r ) where 1 < ij < d, g 1 denotes the monomial g 1 = ■ ■ • gi r G U® r . We will also 
denote / = (i r , . . . , i\). Let 1Z C U® 2 be the set of quadratic relations. Choosing a basis of U® 2 /1Z 
is equivalent to choosing a set C [1, d] 2 such that 

9i9j = E m i,j k ' l 9kgi {i,j)eS( 2 ), 

(fe,0e5( 2 ) 

and fn i ^' 1 = $ij k ' 1 whenever G S^ 2 \ Further, consider the multiindices 

S (r) := {(h, ...,i r ): (ij, i j+ i) € S (2) , 1 < j < r] . 
Finally, let ir be a representation of the algebra on C fc , and we denote Vi = ir(gi) for all i. 

Lemma 2.3. Assume that for all r > 3, the set {v 1 : I G S^} is linearly independent in Mk- 
Then the following intersection property holds: 

1-2 

(2.6) Qi = P| (c d )® x ®g 2 ® (c d )®('- 2 - x ) 

x=0 

for all I > 2. 

Proof. By induction, it suffices to prove 

Gr+1 =Gr®<C d nC d <g>g r =: Gr 

for r > 2. 

Part I. G r +i Q Gr- 

A vector (j) belongs to G r +i if there is a matrix B for which 0j = Tt^Bv 1 ) for / G [1, d\ r+1 and 
it is in G r if there are (Rk)k =0 and (£;)™ =0 with 0/ = Tx(R il v i2 ■ ■ ■ v ir+1 ) = Tr^ • • • v ir L ir+1 ). The 
inclusion C follows by choosing = Bv^ and Li r+1 = Vi r+1 B. 

Part II. G r +i 5 

First, we note that the linear independence condition implies that {g 1 : I G S^} is linearly 
independent at the level of the algebra. Indeed, if it were not the case, there would exist a non-trivial 
relation ajg 1 = 0. In the representation, this implies that = tt(J2i a i9 L ) = J2i a i" 7r (</) = 
^jajv 1 which is a contradiction. Furthermore, spanjg^ : / G S^} = V® r , as every monomial 
g 1 with (ii,ii+±) G S^ 2 ) can be expressed as a linear combination of monomials with (juji+i) G 
S^ 2 \ This, and the linear independence, imply the existence of a unique matrix M^, of size 
d r+1 x |S( r+1 )|, such that 

for all /, where M 1 j = 5 1 j whenever / G S^ r+1 \ 

Now, (j) G Gr +1 iff E/e[i,<+i HTriBv 1 ) = for all 5, i.e. 

/e[i,d]'-+ 1 
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On the other hand, if <f> 6 Q r , then also 

(2.7) <f>j= Ml J^J- 

Hence, 



E ^ v= E E M Vx E M V= E {M*M ( j>)jv j , 

ie[i,d] r+1 ie[i,d] r+1 KeSir+V Jes( r +v Jes ( ^ +1) 

so that Yli < Pl v = implies (M*M)cj) = 0, by the assumption. Since the matrix M has the 
following block structure: 

'If 



M ->N 

it follows that M*M = 1 + N*N > 1, so that det M*M / 0. In particular, (M*M)<p = implies 
cj)j = for J G S'( r+1 ), and by (2.7) = 0. Hence, Q^r + iC\Q r = {0}. This, combined with <5 r +i C Q ri 
implies that G r +i = Q r . □ 

The last general result we recall here is about a uniform lower bound for the spectral gap, see |13j . 
Assume that Hum > is a positive Hamiltonian corresponding to a translation invariant nearest 
neighbor interaction and such that Ker(H[± ^) = G[\,n] f° r a H N > 2. Let Gn be the orthogonal 
projection onto Gn- If there exists k, an e& < \j\fk and iV = iV(fc) sufficiently large such that 

( 2 -8) gk,N ■= ||C[A^-fc+2,7V+l] (G[l,N] - G[i^n+1]) || < £fe • 

Then the spectral gap 7^ of Hnm is bounded below, 

(2.9) 77v > (l-e^) 2 , 

uniformly on the length of the chain, where 7^ is the spectral gap of the Hamiltonian on the chain 
of length k. 

2.2. The MPS representation of the PVBS models and gap estimates. We now turn to 
the analysis of the PVBS systems using their matrix product representation. In this case, d = n+1. 
The quadratic relations satisfied by the generating matrices {wj}f =0 are given by 

(2.10) = e^XiX^VjVi , i ^ j , 

(2.11) vf = 0, ijLO, 

and Ao = 1 corresponds to a choice of normalization. Recall that a vector <f> G G2 iff Yli j=o 4>ij v j v i = 
0, where 4>ij = (e^ (g) ej, 0). Comparing ( 2. 10|2. 11 ) with ( 2.l|2.2 ), it is clear that the vectors ^ and 
4> n span G2 if v^j ^ for all pairs i 7^ j. 

Lemma 2.4. There exists a 2™ -dimensional representation of the commutation relations (2.1(^2A1 ) 
with v^j 7^ for < i < j < n and Vq 7^ 0. 



Proof. In an arbitrary basis of C 2 , let 



a 



















1 



and 

'1 0\ „ f e w 'j/ 2 



Xi) p%3 \ 1 
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The matrices given by 

n 

(2.12) «0 = ®P^» 

i=l 

i—l n 

(2.13) Vi = (g) ® a+ ® (g) P ife 4 , 1 < % < n . 

j=l k=i+l 



are of dimension 2 n and satisfy the commutation relations. Indeed (2.10) foll ows fr om a + di = 
\idicr + and PijO~ + = e l9i ^ 2 a + P^ as well as [Pjk,di] = for all i,j,k. Eq. (2.11) is a direct 
consequence of (cr + ) 2 = 0. □ 

Lemma 2.5. For a,b G Z, 6 — a + 1 = iV > 2 ; £/ie ground state space of H[ a ,b] * s eguaZ to £/at. 

Proof. The intersection property for the spaces £/at follows from Lemma |2,3[ applied to the PVBS 



algebra (2.10, 2.11). In the notation thereof, we choose 

5( 2 ) = {(0, 0)} U {(k, I) : < k < I < n} , 
with m M fci j = for all (k,l) G S (2) if i / 0, and = o'-'^c^"- A,A ; ; . Further, 

5( 3 ) = {(0, 0, 0)} U {(0,0, k) : k = 1, . . . ,n} U {(ijk) : < i < j < k < n} . 

Note that for r large enough, there is a bijection between and the set of subsets of {1, ... ,n}, 
so that 

dim{v J : J G S (r) } < 2" 

for all r and it is equal to 2 n for r > n. The linear independence of this set follows immediately 
from the tensor product structure of the representation and the independence of any di with a + . 
Now, by definition of the interaction and the intersection property, we have that 

Ker(# [a6] ) = G N 

for N > 2. i?r a w being the sum of projections, it is a non negative operator so that C/jv is its ground 
state space. □ 

We emphasize the fact that the intersection property refers only to the spaces Qn themselves. 
The matrix product representation of the vectors therein is only a convenient tool to describe them. 
As Qn = Ran(rAr), the ground state spaces can be explicitly described. It is easy to see that any 
product Vi N ■■•Vi 1 in which an index j / appears more than once vanishes. Indeed, using the 
commutation relations, it is possible to bring the two Vj's next to each other, and recall that t>| = 0. 
Moreover, for any subset S = {si, . . . , s m } C {1, . . . , n}, there exists a vector ip^ a b ^ which is a linear 
combination of all basis vectors of the form 

e ® • • • ® e Sl ® • • • ® e Sm ® • • • ® e . 

where each e Sl , . . . , e Sm appears exactly once. Explicitly, 

x,6] 



(2-14) rf aM =T N (B s ), 



where 

(2.15) £ 5 = P®^- 1 ) 0o--® p®( s 2-si-i) ® a" ® . . . ® o-- ® p^" s ™) 
and P = a + a~ . For example, 

b x-a+l 

(2.16) = Yl ( eWMX i) X a e ® • • • ® e ® ej ® e ® • • • ® e 
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where e% is at site x in each term of the sum. Note that ip^ a y is unique, given the subset S. 
Indeed, the commutation relations determine all components of the vector up to normalization. 
More generally, 

< a E) ■ ■ ■ e Sl 



e 0,^[a,b] 



Finally, the vector corresponding to m = is the simple product 



n 



it 



2.r , 



e ■ • ■ e • 

As S ^ S' => (V'u &] ' V'ja 6] ) = 0) the dimension of the span of such vectors is |"P({1, • • • ,n})\ = 2™ 
for all N > n. Since this is also the dimension of {v J : J E S^ N '}, we have that 



span{^:5GP({l,...,n})} 



/AT 



Since the initial space of Tn has dimension 4 n , its kernel is non trivial for all N. In the language 
of matrix product states, the fact that Tn is not injective for all chain lengths implies that the 
representation is non minimal. And indeed, the minimal representation here is one dimensional, as 
the state on the infinite chain is a product. Here, the larger algebra allows for a richer structure of 
the state when boundaries are present. 

The natural interpretation of the ground states in finite volume is that of a product vacuum ft 
upon which n types of particles can be added, but at most one of each type. Each particle is bound 
to one of the edges, to the left edge if Aj < and to the right edge if Aj > 0, as can be read from 
the exponential decay of the wavefunction's amplitude as the particle is placed farther away from 



its boundary, see (2.16). 



Recall that E is the transfer operator on A4k associated to the matrices vi generating the PVBS. 
Lemma 2.6. For any choice {Aj}™ =1 such that for all i, the spectrum o/E is given by 



Spec(E) 



II A^e^ 1 ^ : n 3 



0,1,2 



l<j<n 

Moreover, the largest eigenvalue Yli-\ >i Af is non degenerate. 

Proof. The proof is by induction on the number of particles. If only one particle is present, say 
with index 1, the spectrum is easily seen to be exactly {1, e i001 Ai, e^ 01 Ai, Af }. Let now d be 
a non empty subset of {1, ...,«}, let j £ S, and denote = S \ {j}. An arbitrary matrix 
M E At 2 |S| can be decomposed with respect to the tensor factor corresponding to j as M = 
M + a + + M_ a~ + M P P + M Q Q, where P = a + a~ and Q = 1 - P. A direct calculation 
yields 



E S (M) 



-W 0j 



•i AjE 5J (M+) ® cr+ + e iy <y \jW J (M_ ) 



W i 



a 



+ E SJ (M P )®P + 



'v^yMp(vl J ) + X]E S3 (M Q ) Q 



where K b = ^2 ie s( v iT ' v i an< ^ v i are the matrices of (2.12 



2.13), truncated to the \S\ factors of 



S. The spectrum of E 5 therefore contains the spectrum of E SJ multiplied by e ±l6,0j Xj . Moreover, 
the eigenvalue equation in the P, Q block can be cast as 



(2.17) 



E SJ - /j, 



0, 
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the determinant of which is simply det(E s ' : ' — /i) det(A|E s ''' — /i). Hence, fj. G Spec(E ,s ) if either 
\x G Spec(E 5J ), or fi/X'j G Spec(E s ' : '). In summary, we have that 



(2.18) 



Spec(E' b ') = {ly^A^e-^A^} • Spec(E 



i0 oj 



which closes the induction. 

Clearly, the largest eigenvalue is n« a^>i ^i- ^ s non-degeneracy follows from the recursion (2.18) 
and the fact that Xf ^ 1 for all i. □ 



In particular, E A is injective for all S. It follows from (2.17) that Mp = 
factor of the corresponding eigenvector. Hence, if there is an index i with A, 



whenever A| is a 
> 1, then the right 

eigenvector for the largest eigenvalue is not strictly positive. Similarly for the corresponding left 
eigenvector whenever there is a j for which A,- < 1, and it is therefore not faithful as a state on 
Mk- In both cases, the representation is not minimal. 

Rescaling E by the largest eigenvalue of the transfer operator yields a completely positive map 



satisfying the conditions of Section 2.1 and generates a unique translation invariant state on the 



chain algebra. We shall now describe the limiting state as well as the ground state spaces of the 
half-infinite chains. First, we consider the asymptotics of the norms of ip£ *n as N — > oo. It is 



straightforward to observe that the same holds for ip^ 
\ % 1. 



-JV+1,0] 



N] 

by replacing the conditions Aj 5s 1 by 



Lemma 2.7. Let ipfj be the vectors given by (2.14)- There exists a finite constant C s > such 
that 



lim _ 



N I 



i£S:\i>l A i 



2N 



where the denominator is equal to 1 if the product is empty. 



Proof. Recalling (2.5), we have that 



(2.19) 



E 

p,q=l 



Xp M [N) [(B s r\ Xp )(x q \B S ) Xq 



J2^Tr{L k (B s )*R k B s ) 

k=l 



where {Xp}p=i is an arbitrary orthonormal basis of C k , and R k , L k are the right and left eigenvectors 
of the transf er ope rator for the eigenvalue t k . Since all Rk's are simple tensor products and with 
B s given in (2.15), it is immediate to observe that (B s )* RkB s = whenever R^ is not diagonal. 
Moreover, B s ^P j B s P j iff j $ S. But PjR k Pj = if \) is a factor of t k . Hence, the only terms of 
the sum that do not vanish are those corresponding to the eigenvalues of the form t k = Y\ - X? where 
the product runs over any subset of S. In the N — > oo limit, the norm is therefore dominated by 
the eigenvalue rije5 A,>i • Note that if {i G S : Aj > 1} = 0, all summands vanish exponentially 



':Ai>l 

but the one corresponding to the eigenvalue 1. 



□ 



Let 



denote the expectation value of a local observable A in the state tjj. 
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Proposition 2.8. Let A be a local observable, namely A G B(7ii a a) for some a < b G 7L. Then, 
for any S G 7>({1, . . . ,n}), 



lim 



/or any x > 0, lim (^4) 
/or any y > 0, lim (.A) 



[a — a:,6+y] 



[a — x,b+y] 



[a — x,b+y] 



,b+y] 








where Sl = {i G S : Aj < 1} and Sr = {i G 5" : Aj > 1} ; urai/i the convention "0... = 0. 

Proof. For simplicity, we consider the case of S* = {1}, and we drop the index 1 everywhere. The 



general case follows similarly using Lemma 2.7 We have 

b+y 

( A )^ la - x , b+y] = \\4>[a-xMv] IP 2 Yl e~ Wii ~ j) X i+j (eo • • • ei ■ ■ ■ eo, Aeo ■ ■ • ei • ■ ■ e ) 

i,j=a— x 

where ei is at position j on the left and j on the right. All summands where one index belongs to 
[a, b] and the other does not vanish. Therefore, 



( A )iP [a - x , b+y] ~ ||^[a-*,6+»]| 



^2 ^ ( e o • • • e , A [ajb] e ■ ■ ■ e ) 

i£[a,b] 



+ e W(i j) A i+j <e ---ei---eo,A M eo---ei---e ) 

i,j€[a,b] 



T!+T 2 



where A[ a b ] is the restriction of A to the interval [a, b]. Now recall that H^a-aj.b+y] || 2 = Y^i=a-x ^ 2 '' 
which diverges as x — > oo if A < 1 and as y — >• oo if A > 1. Since the interval [a, b] is fixed and 
finite, lirn rj2/ _ >00 \T%\ = 0, independently of the value of A, and further 



lim 

x,y— »oo 



(A) 



^[a-:c,i>+y] 



(A) s 



< 



e ---eo,A[ o6 ieo---e ) lim 



1 



|"0[o-a;,6+y] | 



+ lim |T 2 | 

x,y—>oo 



0. 



i^[a,b] 

which proves the first statement of the proposition. We now assume that A < 1, the other case is 
the same after exchanging the roles of the left and the right boundaries. As x — > oo, the calculation 
above continues to hold with y fixed. Finally, if x is fixed, lim^oo h/'[o-a;,6+y] = (1 — -^ 2 *) 
that both liniy^oo |Ti| and lim„_ s>00 \T 2 \ are finite, non-vanishing constants. 



2i\-l 



so 

□ 



Summarizing, there is a unique ground state in the thermodynamic limit, i.e. on Z, which is a 
simple product state fL There is no correlation between spins at different sites of the chain. For each 
of the two half-infinite chains, the ground states are determined by their particle content, with only 
left particles appearing on the right infinite chain with a left boundary and symmetrically for the 
left infinite chain. The binding of the particles to the edges is reflected in the non trivial spectrum 
of the transfer operator. There are 2 nL , resp. 2 nR , ground states on [l,+oo), resp. (— oo,0]. As 
can be read from the proof, the convergence to the limiting ground states is exponential, with rate 
at least 



(2.20) 



e = max {min {Aj, A, X } : 1 < i < n} 
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Note that it would be wrong to conclude that these particles are massless. In fact, there is a 
mass gap for each of them when they are in the bulk, see (2.3) and its proof below, but the binding 
energy to the edge compensates exactly for it as long as no other particle of that type is already 
bound. 

The rest of this section is devoted to the proofs of Proposition 2.1 and Theorem 2.2 We start 
with the lemma below which is a simple but crucial property relating a ground state of the long 
chain to the ground states of the two subsystems obtained by truncating the chain. 

Lemma 2.9. For any S C {1, . . . , n} and < k < N, 



ip N (S) = ^2 T Nik (S k )ipN-k{S\S k ) ®tpk(Sk) , 



s k cs 

\S k \<k 



where 



\TN,k(Sk 



N-k 



Proof. The decomposition follows from the intersection property and the fact that any particle of 
S can appear at most once in ip]y(S). The position of any particle of S k in ip]y(S) is N — k steps 
further from the left boundary than in ^k(Sk), so that each of them is associated with an additional 

□ 



factor Af k , up to phases. 



Proof of Proposition 2. 1 , Part i: Lower bound. First we note that the intersection property implies 
that 

G[1,N] — ^[ljJV+l] = (A ~~ ^[l.iV+l])^!^] , 

so that a vector $ the range of Gum — Gttjv+i] is characterized equivalently by 

$ e G[i,N] n <?[i>+i] • 

Then, 

2 



2 

9k,N 



SUp ||G[7V-fc+2,A r +l] 

$eS [1>N] nG^ N+1] 
11*11=1 



$||" = sup (G[jv_fe + 2,iv+i] 
11*11=1 



$ G £/[i,jv] implies <i> G £7[i,jv-ft+i] by the intersection property. Moreover, P^jv-fc+l] commutes 



with G[7v_fc+2,Af+i] so that 



G 
Hence, 

(2-21) 5-|,jv < 



[AT-fc+2,Ar+l] 



[l,N-k+ 



i] n ^[tv_ 



fc+2,iV+l] 



'[l,JV-fe+l] 



<{N-k+2,N+l] 



± su p m E E 

*ee [1JV+1] , ||*||=1 p,p',5c{l,...,n} i=0 



fe,Jv(P, POMS' 



W_ fc+1 (P)||[|^ fc (P')[|||^(5')|| 



(^ N - k+1 (p)®MP'),^N(S) 



where 4>k,N(P> P') = vpN-k+i(P) ® ^(P'), <&) and tpj^(S,i) = (iPn(S) (g)ej,$). The normaliza- 
tion ||<3?|| = 1 and the orthogonality of the vectors {^(SOIs implies that |( / 97v(5', i)| < 1 and 
l^ivipP')! < 1. Moreover, $ G ^ )iV+ i]' so that 

<p N (S',i) 



(2.22) 



5' i 



\\Hn+i(S) || 
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for all S. By Lemma |2.9[ this reads 

ip N (S',i) 



S',i jeSu{o} 

<PN(S 3 ,j) 

jesu{0} 



N (S')\\U N+1 (S)\\ 

T N+ i A (j) . 



Tn+i,iU) (i>N(S j ) <g> ej,ip N (S') <g> e*) 



nw(^)ir 



Denoting by S L/R := S n Au R and observing that \\ip N (Si)\\/\\ip N+1 (S)\\T N+1)1 (j) = O(Af) if 
i G by Lemmas 2.7 and 2.9 we obtain 



> ip N {b J ;j)— — — T N+ltl (j) 



<Ce N L 



(2.23) 

ie5 fl u{0} 

for any S (it vanishes if Sl = 0), where C is a constant and 

el = max{Aj : i G A^} . 

Furthermore, 

(ip N -k+i(P) ®Tpk(P'),^N(S) ®ei) = lieP'^s=PuP' lT kA i ) T N,k-i(P n )\\^N-k+i(P) 
Hence, 



yfc— i 



(F 



/«M|2 



5fc,7V < sup ^ ^ 



P,P' i&P' 



t> k , N (P,P')<p N (PUP"',i) 



^N-k+l(P)\\Hk-l(P 



■N 



(P U P n 



-T fe) i(£)Tjv,fc-i(-P' 



Asymptotically, for iV > 2k — 1, the coefficients are 



lljeP'nA L A i 



and hence vanish as C±e with e defined in (2.20) in all cases but P n Ar 
particular i G AjjU {0}. The latter terms can be recast as follows 

qIn < Cie k 



P' n A L , and in 



+ C 2 sup 

* s 



ies B u{o} 



-TkAi] 



i> J v- fc+1 (5\5 fl )||||V> fc - 1 (5j ? )|||r jiV|fc _ 1 (5j ? )| 
ll^v(#)|| 



It remains to observe that the second quotient is uniformly bounded and that the sum in the 

mal 



absolute value is exponentially small for any $ by the orthogonality condition (2.23), as 

Un(S' 



bk - l{S ^-T k ^(l + 0(e k ^)) 



for any N > 2k — 1. Finally, since £l < s, there exists a finite constant K such that 

g 2 k>N < Ke k - X < l/k . 

for k sufficiently large and all A" > 2k — 1. This yields a uniform lower bound on the spectral gap 



by the remarks in Section 2.1 



Part ii: Upper bound. We shall use the variational principle 
7 (A,G,A0= inf (H lab] )= inf (H lab] ) . n = inf wM^r^tluo ■ 

First, we observe that subspaces with a fixed number of particles are orthogonal to each other and 
invariant for the Hamiltonian. Let us consider H^ a ^ , the Hamiltonians restricted to the subspace 
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containing exactly one particle of type i and none other. In order to study its spectrum, we center 
and rescale it to 



K 



a,b] - = 



-(Ai + A^Xi^-l) 



which has matrix elements 



'e-^Sa+ifl + e Wi0 6 x - liy a < x < b , 
A" 1 ^,?, + e~ ie,0 5 a+ i,y x = a , 
e i6i0 6 b -i ty + XiS bty x = b. 



On a finite chain [a, 6] of length N, let w be a variational vector of the form w(x) = exp(i/cx) for 
a < x <b and any k € M. Then, 



1 10, tff^ = (AT 1 + At) + (iV - l)£(fc) 
where £?(&) = 2cos(#io — fe)< Moreover, 



a,b] 



W 



\G[a,b] W \ 



N 



[a,b] I 



=: iV - CW(*0 



where & , is the ground state given in (2.16). In terms of the actual Hamiltonian, we have 

[ w > H Ub] w ] 



(2.24) 



N — 1 
N-C N (k) 



E(k) 



Since both | (-0[ a 6 i , w ) | and 



are geometric series, of ratio Ajexpi(/c — 6io) respectively A?, 
the quotient Cjv(/s) is convergent as iV — )• oo for any Aj 7^ 1. Since (2.24) holds for any the least 
upper bound is obtained for k — #io = 0, in which case E(6io) 



c^o) = ^^2< (1+At) 



(1-Ai)(l + A 



AM 



for all TV. Hence, 



1 - G\ 



< 1 



a, 6] 



10 



A.7 1 + A, 



11 — A 



1 + 



2 and 

C; 



C, 



1 



N-d 



□ 



With the normalization Ao = 1, the assumption that Aj 7^ 1 for alH 7^ is essential in the proof 
of the existence of a spectral gap. In fact, the variational upper bound shows that the gap closes 
indeed as |Aj — 1| — > for some i 7^ 0. 

The upper bound of Proposition (2.1) is sharp in each invariant subspace, in the sense that 
1 — 2/ (A" 1 + Aj) also a strict lower bound for finite N, and therefore gives the exact value of 
the gap in the thermodynamic limit for the Hamiltonian restricted to the one-particle spaces. We 
continue working with K % for which the spectral picture there is inverted: The minimal eigenvalue 
for the Hamiltonian H l is the maximal one Aj + A" 1 for its transformed version K l , and the 
gap bound 1 — 2/(Aj + A" 1 ) for H l corresponds to 2 for K l . Let us write the eigenvalue equation 
K l N+2 ip = Eijj using the transfer matrix 



T(E) = 

The boundary conditions 



1 







e i0iO 

1 



E -1\ fl 
1 Ho e Wi0 



^(l)=e ie ^(E-\7^(0), VW 



-iSiO 



(E-Xi)^(N+l), 
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yield the equation 

(2.25, •"""""(f oTCV"" 1 * 

on the chain of length N + 2. Since 



1 ; " \U N -!{E/2) -U N „ 2 {E/2)J 
where Ujy(x) are the Chebyshev polynomials of the second kind, ( |2.25 ) reduces to 

(E - Aj) [C/iv(^/2) (£ - Ar 1 ) - U N ^(E/2)] = U N ^(E/2) (E - Ar 1 ) - U N . 2 (E/2) , 
and further to 



(2.26) U N {E/2) 



+ 



(E-Xi) {E-X; 1 ) 



U N ME/2)- {E _ x _ 1){E _ Xi) U N ME/2). 



If E = X{ + A^ 1 > 2, this equation is precisely the defining recursion relation for the Chebyshev 
polynomials, so that Aj + A" 1 is an eigenvalue for all N, and indeed the largest one. Moreover, as 
Un(1) = N + 1, we also note that E = 2 does not belong to the spectrum for any finite N, as (2.26) 
reads 

l _ 3-jXj + Xr 1 ) 
5 - 2(Xi + X; 1 ) 

whose only real solution is Aj = 1, which is excluded. Now, we can use the eigenvalues of the 
transfer matrix 

t ± = ±(E±VE*-£ 

to diagonalize it away from \E\ = 2 and obtain the following form of (2.25): 

(1 -t.{E- Xi)) i;^ = eW+V^tZ {(E - Ar 1 ) - t _)) An 

(-1 +t+(E- Xi)) Am = eW+Wni* (-(E - XT 1 ) + 1+)) Ar, 

If E / Xi + A" 1 , the system yields 

l-UE-\ i )t+-(E-\r 1 ) = /t+^ A 
l- t+ (E- Xi)t- -(E- Xr 1 ) \t- 

The left hand side is strictly decreasing for EG [2, Xi + A" 1 ) and its limiting value at E = 2 is 1, 
so that 

l-t-(E-\i)t + -(E-\ri) = /t+N v 

l-t+(£7-Ai) AT 1 ) V*- 

which is a contradiction as £+ > t_ for E > 2. Hence, there are no eigenvalues of K l N+2 m the 
interval [2, Aj + A^ 1 ). Therefore, if rf(N) is the spectral gap above the ground state of H l N , we 
have that 

inf (f(iV)) > 1 — r . 

W ^ ! > ~ Xl + XT 1 

We further conjecture that 

f 2 

min < 1 r : 1 < i < n 

I Ai + Ar 1 " " 

is in fact the exact gap of the full system in the thermodynamic limit. 
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Proof of Theorem 2.2. We recall that if two models are in the same gapped ground state phase, 
then for any lattice under consideration their ground state spaces are mapped onto each other by a 
quasi- local automorphism. Suppose that n L ^ n 2 L . Then the dimension of the ground state spaces 
of the two PVBS models on the half-infinite chain with a left boundary are not equal, and the two 
sets cannot be related by an automorphism. The same holds for n l R ^ n 2 R and an infinite chain with 
right boundary. Hence, the condition (n^,n R ) = (n L ,n R ) is necessary for equivalence. In order to 
show that it is also sufficient, we construct a smooth interpolating path of gapped Hamiltonians and 
conclude by [I]. The parameters of the two PVBS models are denoted by ({Aj}, {%}), respectively 
({/Uj}, {tpij})- By relabeling the basis elements ej, we can assume that 

Al < A2 < • • • < Ki L < 1 < A„ i+ i < • • • < X nL+nR . 

Let (i 1, . . . , i n ) be the sequence of indices for such that 

Mil < • • • < IHn L < 1 < Mi„ L + l < • • • < Vin L +n R • 

Let u be the unitary map defined by 

uei. =Bj, ue = e , 

and u(s) 6 SU{n + 1) a smooth path for s £ [0, 1], such that u(l) = 1 and u(0) = u. Finally, let 
U(s) = <3) x=a u(s). I n the same interval, let 

Xj(s) = s\j + (1 - s)/ij 3 

0jk(s) = s6 jk + (1 - s)(fi. ik 

and h(s) be the PVBS interaction constructed from these s-dependent parameters. Then the 
Hamiltonian 



H [aM (s) :=U(s) 



Em 



S)x,x+1 



u(sT 



satisfies all the properties required for the construction of a quasi-local automorphism. Clearly. 

[aj] aIld %,&](!) = H [al] 

IS 



H[ a ,b](Q) = Hfc'b] and H[a,b](l) = ^tfb}- Moreover, the linear interpolation of the parameters 



generates a smooth family of normalized ground state vectors ip^ a 6 j(s), which remain orthogonal to 
each other for s G [0,1]. This implies that the path h(s) is smooth as well. Since all paths Xj(s) 
interpolate linearly between parameters that are either both larger or both smaller than 1, we have 



that \Xj(s) — 1| > for all j and s £ [0, 1]. Hence, Proposition 2.1 ensures that the Hamiltonian 
Ylx~^a h{ s )x,x+i is gapped, uniformly in the length \b — a + l\. Finally, the unitary conjugation 
preserves the spectrum and therefore the gap, and since U(s) is local, the Hamiltonian H< a u(s) has 
a nearest neighbor interaction. □ 

3. The ground state phase of the AKLT model 

In the previous section, we introduced a family of simple models that can be classified according 
to their boundary behavior. They are all gapped with a unique ground state in the bulk, but they 
support a number of additional edge modes in their ground state spaces as soon as boundaries are 
present. In the bulk, these states become gapped excitations. 

We believe that these characteristics are the key for the classification of gapped phases beyond 
symmetry breaking in one dimension. We now illustrate the use of the PVBS classes with a physi- 
cally relevant example: The AKLT model [12] . It is a one-dimensional spin-1 chain, namely T-L x = C 3 
carrying the three dimensional representation of su(2), with generators denoted (S , S 2 , S 3 ). Its 
Hamiltonian has a nearest-neighbor interaction 



H AKLT = E 



1 1 2 1 

- (S x ■ S x+ i) + - (S x ■ S x+ i) + - 
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corresponding to the projection onto the spin-2 subspace of every pair of neighbouring spins C 3 (S>C 3 . 
As such, Haklt > 0. Its unique ground state in the thermodynamic limit has a matrix product 
representation identified in [18] described by the set of matrices 

in a fixed basis {xi)i=± of C 2 . They satisfy the following commutation relations 

(3.2) WqW- = —W-Wq , wqw + = —w + wq , w± = , 

(3.3) W-W + + w + W- = —2wq . 

For notational clarity, we denote here the transfer operator and the matrix product map on M2 by 
T and T^v- The uniqueness of the thermodynamic ground state follows from the fact that 1 is an 
isolated eigenvalue. Moreover, and unlike the PVBS situation, there are short range correlations 
between sites in the thermodynamic limit corresponding to the triply degenerate eigenvalue —1/3. 

The four-dimensional ground state space of a finite chain, namely the range of T^r, has the 
following picture: It is isomorphic to the tensor product of two Bloch spheres that can be thought 
of as sitting at the two boundaries of the chain. Despite the superficial simplicity of this description, 
the ground states have a complicated structure that can be understood as follows. In the product 
basis, the only words with non vanishing coefficients in any ground state are those with alternating 
'+' and '— ', separated by an arbitrary number of '0'. The four ground states differ by whether 
the first non zero letter is a '+' or a '— ' and whether they come in equal numbers or with one 
additional '+' or '— '. In particular, they do not differ in the bulk and they all converge to the 
same weak-* limit point in the thermodynamic limit. Moreover, the system on the half-infinite 
line has a two-dimensional ground state space isomorphic to a single Bloch sphere. For a more 
detailed description of the model, but in the slightly different language of valence bond states, we 
refer to the original article [12] ■ Furthermore, \W\ I2U] identify a 'hidden string order' as a rigorous 
description of the 'dilute Neel order' heuristically introduced above. 

With this quick introduction, we can now state the main result of this section. 

Theorem 3.1. The AKLT model belongs to the PVBS phase with til = nR = 1- 

In other words, there exists a quasi-local automorphism of the observable algebra that maps the 
ground state spaces of the AKLT model on finite, half-infinite and infinite chains onto those of the 
PVBS class with one left and one right boundary particle. In particular, the intricate bulk state of 
the AKLT chain is equivalent to the simple product state common to all PVBS models. 

The importance of the locality property of the automorphism comes into sharp focus when one 
considers the ground states on the infinite and half-infinite chains. In finite volume, any unitary will 
preserve the dimension of the ground state space. The non-local unitary transformation introduced 
by Kennedy and Tasaki in |20] to reveal the hidden string order transforms the 4 AKLT ground 
states into 4 translation invariant product states, which leads to 4 distinct bulk ground states in the 
thermodynamic limit. This example shows that non-local unitary transformations do not preserve 
the structure of the bulk ground state(s) and clearly would not be useful to classify gapped ground 
state phases. 



To prove the quasi-local automorphic equivalence of Theorem 3.1 it suffices again to construct 
a smooth path of uniformly gapped Hamiltonians interpolating between the two models. We shall 
first construct a path of interactions, show that the Hamiltonians have matrix product ground 
states and use this property to prove the existence of a gap. 

In |llj . we presented the result without complete proof. In order to make the comparison with 
the current article easier, we note the change of notation from A(s) and /i(s) there to f(s) and g(s) 
here. 
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3.1. A path of Hamiltonians. In order to obtain a smooth path of Hamiltonians, we define a 
smooth path of nearest neighbor, translation invariant interactions interpolating between the AKLT 
interaction and the PVBS interaction for i%l = Ur = 1. Both extreme points being projections, 
it suffices to define a smooth path of vectors, and to use their span as the range of a path of 
projections. Let so be defined by sin(so) = y / 2/3. For s G [0,so], let h(s) be the projection onto 
the five dimensional space spanned by the vectors 

e+ ® e + , e_ ® e_ , e+ ® e + /(s)e ® e+ , e ® e_ + /(s)e_ ® e , 

( 3 - 4 ) , v sin(s) , 2 , . 

e+ ® e_ + gr(s) — ^^eo ® e + J (s)e_ ® e+ , 
cos z {s) 

where 

(3.5) | 5 (s)| 2 + |/( S )| 2 cos 2 (s) = l. 



Lemma 3.2. Assume that s i— )• /(s) is C 1 ([0,so]), i/iai |/| < 1 and let g(s) be a solution of (3.5). 
Then, for any a < b G Z, i/ie map 

6-1 

s 1 — > H [aM (s) = ^ K s )x,x+l 
x=a 

is of class ^([O^o];^^ 6 -^ 1 ))). 

Proof. If /(s) is C 1 ([0,so]) then all of the vectors above are strongly continuously differentiable 
in [0, so], and since they are orthogonal to each other for all s, the projector h(s) is continuously 
differentiable in [0, so] in the strong operator topology. Hence s h-» H^ a b ^(s) is strongly C 1 , and 
therefore also in norm for all a < b G Z. □ 

Let now / be an monotonely increasing real function satisfying the boundary conditions 

/(0) = 1/V2, /(s ) = l, 



At s = sq, the interaction is the AKLT interaction. At s = on the other hand, the range (3.4) of 



the interaction corresponds precisely to the vectors ( 2.1|2.2 ) orthogonal to the ground state space 
of the PVBS model with parameters 

A + = v / 2, A_ = l/V2, 

and 

= e 0+ = 9 - = it. 

Hence, the Hamiltonian at s = is a PVBS model with i%i = ur = 1 and the path of Hamiltonians 



defined in Lemma 3.2 interpolates as announced. It remains to show that these operators have a 



uniform spectral gap above the ground states. 

3.2. A path of matrix product states. In this section, we show that the Hamiltonians can be 
constructed pointwise for s > as the parent Hamiltonians (see [13]) of an interpolating family of 
matrix product states. For s G [0, sq], we define the 2x2 matrices 

(3-6) w-(s) = (° , w (s) = -cos(s) (J _ / ° (s) ) , W+ (s) = (j^ J 

where /(s) and g(s) were introduced above. For all s, let T^(s; B) be the completely positive map 
corresponding to these matrices. 



Lemma 3.3. Let f be as in Lemma 3.2 and g be a solution of (3.5). Then for s G (0, so], there exists 
an ergodic purely generated MPS on the spin-1 chain corresponding to the matrices {iOi(s)}i=±,o- 
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Proof. First, we note that Condition (3.5) ensures that W(s) is an isometry, or equivalently that 1 
is an eigenvector of T(s) for the eigenvalue 1, as 

f(s;i) = Wi(sYMs) = 1- 

i=±,0 

In fact, all right eigenvectors are given by 

Ri(s) = l, R 2 ( s ) = a + , R 3 (s) = a~, R 4 (s) = q 2 (s)P - gi (s)Q 

where Q2(s) = (sin 2 (s) / g 2 (s)) Qi(s). The associated eigenvalues are given by ti(s) = 1, t2(s) = 
£3(5) = — /(s)cos 2 (s) and t 4 (s) = cos 2 (s) — g{s) 2 = /(s) 2 cos 2 (s) — sin 2 (s). The assumptions on 
f{s) ensure that 1 is non degenerate for all s 6 [0,so]- The corresponding unique left eigenvector 
is given by g{s) = Q\{s)P + Q2{s)Q. With the normalization, we get 

/ \ 5 f ( s ) 2 / \ sin 2 (s) , , . 

8i(s) = 7T^— t^tt , Q2{s)= . 2 =l-gi(a), 

g(s) + sin (s) g{s) + sm (sj 

and the matrix £> defines a faithful state on the auxiliary algebra M2 for s > 0. Hence, the triple 
(T, 11, q) generates a translation invariant state on the chain algebra, which is moreover extremal 
since the peripheral spectrum of T(s) is trivial for all s G [0, so]. □ 

In fact, the transfer operator T(s) can be completely diagonalized 

4 

f(s; A) = Y^ U(s) Tr(Li(8)A)Ri(s) , 

i=l 

with the left eigenvectors Li(s) given by 

L 1 (s) = g(s), L 2 (s) = R 3 (s) L 3 (s) = R 2 (s) , L A {s) = P - Q = a* . 

The normalization has been chosen so that Tr(Lji?j) = 1. 

The matrices {wi(s)}f =0 form a representation of the following deformation of the algebra (3.2) 
of the AKLT model into that of the PVBS model: 

(3.7) f(s)w (s)w-(s) = -w-(s)w (s) , w (s)w + (s) = -f(s)w + (s)w (s) , w±(s) 2 = , 

(3.8) w^(s)w + (s) + f(s) 2 w + (s)w^ S ) = _ g( g ) sin ( g ) Wo(s) 2 

cos z (s) 



Although the representation ( |3.6[ ) holds for all s, including s = 0, it fails to be faithful at the 
boundary point as w+(0) = 0. A non-trivial four-dimensional representation of the algebra at 
that point was given in Section [2] Interestingly, the following lemma shows that the discontinuity 



in the dimension of the minimal non-trivial representation of the smooth deformation ( 3.7|3.8 ) is 
necessary. 



Lemma 3.4. Let v~, vo,v+ be a two-dimensional representation of the commutation relations (2.1(^2A1 ) 
with n = 2 and A_ 7^ A+. Then v- = or v + = 0. 

Proof. Assume that V- 7^ and v+ 7^ 0. The condition A- 7^ A+ together with the commutation 
relations for ( — , 0) and (+, 0) imply that V- 7^ v+. In two dimensions, any nonzero nilpotent matrix 
must be unitarily equivalent to a + . Hence, we let V- = acr + , v+ = f3U*v-U, with U 7^ 11 and 
(3 7^ 0. Taking the trace of the relation for (— ,+) and using cyclicity yields Tr(u_u + ) = 0. This 
implies (w+)2i = 0, so that U12U21 = 0. Hence, by unitarity, U = 1, which is a contradiction. □ 

At this point, one may also wonder why the commutation of w-(s) with wo(s) and that of WJ r (s) 
with wo(s) involve parameters that are inverses of each other. Indeed, at the PVBS point, A+ and 
A_ can be chosen independently of each other. Here again, it turns out that it is in fact dictated 
by the form of the commutation relations themselves. 
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Lemma 3.5. We consider the following quadratic relations: 

W-Wq = C-qWqW- , W + Wq = C + qWoW + , w± = , 

W-W + = C |_u) + ii;_ + Dwq , 

where all coefficients Cij and D are non zero. If there exists a faithful representation such that 
Wq / 0, then C + o = l/C_o. In particular, 

[w + w- , Wi] = , i = 0, ± . 

Proof. On the one hand, 

W-W + Wq = C-qC + qWqW-W + = C-qC + qDWq + C-qC+qC ^WqW + W- , 

while on the other hand, 

W-W + Wq = C \-W + W-Wq + DWq = C ^C-qC + qWqW + W- + DWq , 

so that C-qC+qDwq = Dwq which yields the claim. □ 
In the PBVS case, D = and the lemma does not apply. 

The condition C +0 = 1/C-o is equivalent to the linear independence of {w 1 : I G S^} at the 
level of the abstract algebra, and ensures that the intersection property holds. 

Lemma 3.6. For < s < so, let 1Zn = Ran(Tjv). For N > 2, the spaces TZn(s) satisfy the 
intersection property (2.6). Moreover, &\m(lZN{s)) = 4. 



Proof. The intersection property follows from Lemma 2.3 In the notation thereof, we choose 
5( 2 ) = {(00), (0+), (0-), (+-)}, so that 

5W = {(0---0),(0---0+),(0---0-),(0---0 + -)} 

and the linear independence of the corresponding {w 1 : I G S^} is immediate from the explicit 
representation of the matrices (3.6). That dim(7?.Ar(s)) = 4 follows from = 4 and again the 

linear independence. □ 

It will be useful to have an explicit description of a well-behaved basis of Gn( s )- F° r 2 < N < oo 
and s G (0, so] let 

(3.9) (»( S ):=r N (s;A%(s)), 
where 

A 1 N (s) = P + q N ( S )Q, A 2 N (s) = a~, A%(s) = , A%(s) = [(-f(s)) N P - Q] , 

with 

t2(s) N -(-f(s)) N (g 1 (s) + Q 2 (s)t 4 (s) N ) 



(3.10) g N (s) 



Now, (2.19) reads 



(3.H) <^( S ),CAV( S )> = E^( s ) 7VTV ( L i( s X( s )^( s )^( s )) 
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for any 1 < /j,, v < 4. With Rj(s), Lj(s) given above, this vanishes for all pairs /i < v, the case 



= 1, v = 4 depending on the choice (3.10) for <?at(s). Moreover when fi = u, (3.11) yields 
(3.12) 



f?i(s) + ^2(^*4(5)^) + 2q N (s)t 2 (s) N + ^(s) 2 (g 2 (s) + £»i(s)t 4 (s 



(3.13) 
(3.14) 



|C^)|| 2 = ||Cto|| 2 = ^i( S ) (l 



9{sr 

sin 2 (s) 



{Q2(s) + Qi( S )U(s) N ) 

+ f(s) 2N U(s) + ^ 2 ( S )i 4 ( S ) iV ) - 2(-f( S )) N t 2 (s] 



N 



These norms remain strictly positive for all s > 0, and the set (C^r(s)) 4 =1 forms an well-defined 
orthogonal set spanning the four-dimensional space Qn(s) for s > 0. In fact, lim^o (?jv(s) exists 
and the four vectors Cn( s ) can ^ e extended by continuity to s = 0. Furthermore, 



lim lim kX(s) 
Ar->oos->0 11 JV 

for all ji. From now on, we define 



lim lim C^rfs) 



lim Qi(s) = 1 
s— >o 



(3.15) 



K N (0) :=span{C^(0+) :// = !,..., 4} 



and in particular, 

K N (0) / Ran(T iV (0)) . 

Indeed, since w + (0) = 0, the range of Tat(0) is two-dimensional. The map Tjv(0) generates the 
ground states of the system on the infinite chain without boundary and the chain with a left 
boundary, but not on the other half-infinite system. With minor modifications, the roles of left and 
right could be exchanged. 

As we already noted in the introduction, this example emphasizes the importance of constructing 
a family of Hamiltonians and ground states, as a smooth family of matrix product maps Tjv(s) does 
not necessarily generate the desired continuity if the parent Hamiltonians are simply constructed 
pointwise. 

Despite being a simple corollary of the above discussion, the following proposition is essential: 
it identifies the MPS constructed from the path of matrices for s > and their limits at s = + 
with the ground states of the Hamiltonians of Lemma |3.2| 



Proposition 3.7. LetlZ]y(s) = Ran(Y7v(s)) for s > andTZ^^O) be defined by (3.15). If H[ ab Us) 
are the Hamiltonians of Lemma\3.S\ with N = b — a + 1, then 



n N ( s ) 

for all s E [0, So] . 

Proof. By construction, and for all s > 0, 



Kerf//, 



0,6] l 



0), 



C2( s ) = -f( s ) 2 9( s ) sin(s)e_ <g> e+ + cos 2 (s)(l + f{s) 4 )eo ® eo — g(s) sin(s)e+ <g) e_ 



(3.16) 



Cl(s) 



-/(s)eo ® e_ + e_ ® eo , 
-f(s)e+ <g> e + e <8> e+ , 
_ <8> e+ - /(s) 2 e+ <g> e_ , 



where we used that q 2 (s) = f(s) 2 . These span the orthogonal complement of the space spanned 
by the five vectors of (3.4), so that 7Z 2 (s) = Ker(/i(s)) = ~Kei{H^ xx+ i^{s)). The statement for all 
N > 2 follows by the intersection property. □ 
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Interestingly, we see that the path chosen here maps the spin-0 singlet ground state of the AKLT 
model to the product state of the PVBS model, while the spin-1 triplet is deformed to the edge 
states. 

The last ingredient needed in the proof of Theorem 3.1 is a uniform lower bound on the spectral 
gap along the path of Hamiltonians. The following lemma, whose proof is very reminiscent of that 
of Proposition 2.1 is the key. Recall that Gn(s) is the orthogonal projection on Qn{s). 

Lemma 3.8. Consider the following family of functions on s E [0, sq] 

9k,N(s) := ||G[Ar-jfc+2,iV+i](.s) (G[i,N]( s ) ~ G^iV+i] ( s )) \\ ■ 

for k eN and N > k. Then: 

i. For all N > 2k — 2 and for all s E [0, sq], there is a C(s) such that 

(3.17) g KN {s) 2 < C(s)e k - 1 

where e = sup sg [ 0jSo ] ( m a x ie{2,3,4} l*i( s )l)- Moreover, s h-> C(s) is continuous on (0,so]. 

ii. Let < / < 1 be fixed. Let < s\, S2 < vr/2 be defined by 

■ 2, , I-/ 2 ■ 2, s / 2 (l~21n/) 

sm = YTp ' sm (fi2) = — TTp — 

and let 5 = minjsi, S2, so/2}. Assume that f(s) = f for s < 5. There exists k\ such that 
k>k\ implies gk,N(s) 2 < 1/k for all s < 5 and N > 2k — 2. 



Proof, i. Pointwise bound. By a direct translation of the steps leading to (2.21), but with s- 
dependent objects and ground state vectors indexed by fj, = 1, . . . , 4, we can write 



(3.i8) 9 „ kW * = jS p ) y, E K ^m^nul(M <cft -*« (s) @ ffW ' c " (s) 8 Cl> 

where 

K(s) := E H [1>N+1] : ||$|| = 1 and G g [1>N] (s) D G^n+I}^} , 

<?W;fc,Jv(s) = (CAr_fc+i( s )®Cfc(s),$) and (p Tji;N (s) = (C^(s)<g)ei, $}. The normalization ||$|| = 1 and 
the orthogonality of the vectors (C^r(s))^ =1 imply that \4>^,i/:k,N(s)\ < 1 as well as |</? T ,i;jV"(s)| < 1. 

If s > 0, scalar products of ground states can be expressed using the transfer operator and its 
spectral decomposition, namely 

4 

<C_ fc+1 ®C£,&®e<)= ^ ^ 1 ^- fc+1 Tr(^A] v LX(^)*^(^- fe+ i)*) , 

p,ij=i 

4 

<C*+i, C]. ® = ^ i^Tr (R p A T N L p wt{A» N+1 f) . 
P =i 

With the second identity, the condition <3? G ^[i^at+i]^) -1 " reads 

( 3 - 19 ) E ^r Tr (^wc^+i)*) = E iaft E Tr • 

r,i \^N\\ T7i \\^N\\ p=2 
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The first identity on the other hand yields 



(Cft- fc+ i®Cg,g) 
HCiv-fc+ill HCjfcll 



1 



lie 



N-k+1 



ma 



E JrJ Tv ^NP^i^kTp^N-k+iT) 



lis. 



P+g>2 
l JV+lii/= 



A" I 



Since the matrices (A^ +1 )f, =1 are linearly independent, there exist coefficients c utJ, a (k, N) such that 



E^(^)(A^ +i y 



(T=l 



Note that c CT (A;, iV) are uniformly bounded in fc, JV. It remains to use this decomposition and the 
orthogonality condition (3.19) in the first line above to obtain 



(3.20) gl N = sup YY, 



1 1Kjv-fc+illllCfe IIIICatI 



Y N) Y $Tr {A N pw*(A N+i y) 

cr p=2 



p+q>2 

Recalling that there exist < C m < Cm < oo such that 

C m <\\(k(s)\\<C M 

uniformly in N and noting that the traces are finite for any fixed s > 0, the fact that the peripheral 
spectrum of the transfer operator is empty implies the exponential decay (3.17) of gk,N(s) 2 for all 
N > 2k — 2. Finally, we note that the case s = corresponds to the PVBS model and has been 
treated independently in the proof of Proposition 2.1 The continuity of C(s) follows from that of 
all elements of (3.20) for s > 0. 

ii. Uniform bound. A priori, the function C(s) has a singularity at s = due to the factors 
sin(s) -1 of the matrices A 3 and A 4 , appearing both in the traces and in the coefficients c u ^ a {k, N). 
If Bk t N(s) is the square bracket in (3.20), we write 

Bk,N{s) = . o- : + r kN (s) 
sin {s) 

by gathering all singular contributions in the first term. The regular part rk n( s ) containing no 
negative power of sin(s), it can immediately be bounded above as 

Kiv(s)| < C T e k - 1 , 

with an s-independent constant C r , for all iV > 2k — 2. For s < s±, and recalling that f(s) is 
constant there, q\ (s) can be expanded in powers of sin(s) 

1 + / 2W-1 



01 W 



1 



J 3=1 V 



sin J (s) , 



and similarly for Q2(s) = 1 — Qi(s). Since the coefficients of a product are given by convolution, 
the expansions of £i(s) 2 , Q2{s) 2 and q\{s)q<2,(s) converge absolutely for s < s\. Hence, the singular 
part of (3.20) has the following convergent representation, 



(3.21) 



Ck,N(s) 

sin 3 (s) 



E Ev 

3=0 \P=2 



s) a 



i + E tpi^w, 

P+g>2 



N-k+li 



sin- 7 3 (s) 
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2 r , 



where the coefficients a,j = a,j(k, N) and bj = bj(k, N) are independent of s and uniformly bounded 
in k, N. The terms j > 3 are regular and summable for s < s% and can be estimated again by a 



fc-i 



constant times e 

Now, we use the binomial expansions of t p (s) n in the < j < 2 terms of (3.21) and reorganize 
the sums to obtain 




E E 



r k-l r AT-fe+l 



sin* 3 (s) 



p+q>2 a, f3: 

j+2a+2/3=i 

k-l,N-k+l 

E E 

p+<j>l a,/3 

j+2a+2/3>3 



r fc-l T JV- fc+ l sin i+2a+2^-3 (s) 

p,a p,p v ' 



Ai = / 2n QjH 1 +r 2 )) i - 

Now, the a priori bound gk,N(s) 2 < 1 implies that Cfc,Ar(s) = Ofc,Ar(sin 3 (s)), and in turn that the 
first line vanishes. In the second line, where all exponents of sin(s) are non-negative, it suffices to 
observe that 



E 

fc=7 



n 



r.fc-7 - ^ U 



< 



P" < n 7 e nN 



Hence, 



JV 

r . sm 



j+2o-3 



< r ((i + /- 2 )iv) 



,r(3-j')/2l e AT(i + /-2) sin 2( s) 



N 

E 

a:j+2a>3 

and similarly for the other terms. As N — > oo, the left hand side decays exponentially for s < S2- 
Summarizing, all contributions to (3.21 ) show a uniform exponential decay in k for N > 2k — 2, so 
that for k large enough, at(s) 2 < 1/k for all s < 5. □ 

Note that the specific assumptions sufficient for point (ii) to hold are of technical nature and we 
have not aimed for the optimal conditions. Any neighborhood of s = in which the gap remains 
open will do, as part (i) ensures the existence of a gap for s > 0. Similarly, the choice of a constant 
/ around the origin simplifies the perturbative calculations in the proof, but is not necessary. A 
general theory of the continuity of the spectral gap for such families of Hamiltonians still needs to 
be developed. 



Proof of Theorem 3.1 Lemma |3.2| introduces a smooth path of Hamiltonians for s G [0, So]- The 
continuity of the map s h-» C(s) on (0, sq] introduced in Lemma Qi) implies that it is bounded 
on [5, so]. Therefore, there exists k$ such that k > k$ implies 

9k, N (s) 2 < C(s)e k - 1 < 1/k s>5, 

for all N > 2k — 2. Together with (ii) of the same lemma, gk,N{s) 2 < 1/k for all k > max{/co, k\} 
and all s 6 [0, so]. Hence, the martingale result (2.9) implies that the path is gapped, 

inf 7 (s,AO>0. 

se[o,s ],AreN 

Hence, by PQ, there exist quasi-local automorphisms a r , for T = Z, (— oo, 0], [1, +oo) on the algebras 
of the infinite chain, of the right infinite chain with a left edge and of the left infinite chain with a 
right edge, and concludes the proof of the theorem. □ 
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4. The SO(2J + 1) models 

As a further example of the classification of gapped phases using the PVBS classes, we consider 
the SO(d) invariant models introduced in [T7], which are one possible family of generalization of 
the AKLT model. The local Hilbert space is d dimensional and carries the fundamental represen- 
tation T> of SO(d). In order to define a nearest neighbor interaction, we consider the irreducible 
decomposition of the tensor product of two fundamental representations into 

(4.1) V®V=l@A@S 

where 1 is the invariant singlet, S contains all other symmetric vectors, and A is the set of the 
antisymmetric vectors. Let P$ be the orthogonal projection onto S, and let the Hamiltonian be 



6-1 



(4-2) H [a>b] = ^(P s )*, x+ i 



By construction it is invariant under the action of SO(d). In the case d = 3, the decomposition (4.1 ) 
coincides with the Clebsch-Gordan decomposition of the tensor product of the tensor product of 
two spin-1 representations of su(2), and P$ is the AKLT interaction. In the sequel, we shall focus 
on the odd case, d = 23 + 1. 

The matrix product representation of the SO(2J + 1) models' ground states can be chosen to 
be the generators {Z a = Z^\a = 0, . . . , 2 J} of the Clifford algebra C 2 j+i satisfying the quadratic 
relations 

(4.3) ZaZp + ZpZa = 28 a pt. 

The fact that the symmetry remains unbroken in the ground states is manifest in the invariance 
of (4.3) under the transformation Z' = ^2 a 0.y a Z a for any O G SO(2 J + 1). The correct normal- 



ization is obtained by setting Z a = jZ a and imposing the isometry condition 1 = ^ Q Z* a Z a . Since 
Z\ = 1 for all a, we have that 

(4.4) | 7 | 2 = ^^. 
V ; 1 2J + 1 

A representation of C2J+1 can conveniently be obtained from a representation of the CAR algebra 
A-(C J ). »4_(C) is generated by 1, <r + and a~ , and the higher dimensional fermionic creation and 
annihilation operators are given by 

(4.5) a* = (-l) Q ® • • • ® (-1) Q ® 1® • • • ® 1 € -4_(C J ) . 



Note that (-1) Q = a z . Then, 



3-1 J-j 



J J 



Q 



(4.6) Z = Y[(2a* aj - 1) = (g)(-i; 

i=i 3=1 

and for 1 < j < J, 

Z 2 j-i = dj + a* , Z 2 j = i(aj — a*) . 

The canonical anticommutation relations of the {apdj}j =1 imply the Clifford relations for the 
{Z a }^-Q. This representation clearly is of dimension 2 J . 

Lemma 4.1. The ground states of (4-2) are MPS generated by the matrices {V^}^, where 



1 

v / 2T+T~ u ' FZ ^ i_ V2J + l^' V2J + 1^' 



Va = ; _ Zo , V23-1 = \ aj , V 2 j = -\l — — - a* 



for 1 < j < J, and the operators a - are given by (4-5) and Zq by (4-6). 
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Proof. We prove that the matrices of the lemma generate the same vectors as { r yZ a }~L . Quite 
generally, the pullback of a change of local basis on the chain tp a ^ 4>/3 = E a ^3a^a to the 
generating matrices v a of a MPS is given by v a h- > vp = Y^ a U/3 a v a - Indeed, 

^2 Tp/3 <8> V/3X = ^2 U P^a ® VpX = ^^a® V a \ 

P a,/3 a 

implies v a = X^(^ T )a/3^/3- The definitions 

Vo = 7-^0 , V23-1 = ajaj V 2j = f3ja* 

correspond to a change of basis given by a matrix U with Uqq = 1 and 2x2 blocks 

1 / ay — iay 
¥f [fa \Pj 

for 1 < j < J. Unitarity imposes that \aj\ 2 = \(3j\ 2 = 2|7| 2 . The claim follows by choosing 
7 = -(2 J + l)" 1 / 2 in accordance with (4.4). □ 



The arbitrary choice of sign is designed to match exactly (3.1) in the case J = 1. With this 
mapping from the Clifford algebra to the CAR algebra, it is now immediate to generalize the path 
of matrix product maps and therefore of Hamiltonians from J = 1 to an arbitrary J, thereby 
introducing a smooth path of Hamiltonians between any SO(2J + 1) model and the PVBS model 
with J particles of both types. If the path is gapped indeed, this places these higher spin models 
into the PVBS classification. 

We start by introducing a deformation of the CAR algebra to relate it to the PVBS algebra with 
n L = n R = J ■ For any complex number A, 

(4.7) a+o-- + \a-o- + = \ Q , 

(4.8) \-a-\ Q = \ Q o- , a + X Q = A • \ Q a + . 

Given a vector A = (Ai, . . . , Aj), we define the twisted creation operators (a*j(A))j =l by 
a*(A) = (-Ai) Q <g> ■ • • ® {-\j-i) Q <g> <J + ® A^ +1 ® • • • ® Aj , 

and 

ao(A) = (-\i) Q ® ■ ■ ■ ® {-\j) Q . 
Clearly, (a*(A)) 2 = = (aj(A)) 2 . From (4.7) and (4.8) we obtain the twisted commutation rela- 
tions, 

(4.9) a*(A) aj (A) + X 2 a j (A)a*(A) = a (A) 2 

(4.10) a*(A)a fc (A) = -A,A fe a fc (A)a*(A) (j ^ k) 

(4.11) a*(A)a (A) = -X jaQ (A)a*(A) 

(4.12) a*(A)4(A) = -A,A^V fc (A)a*(A) 

A higher dimensional analog of the algebraic path (3.7) is easily obtained with the following defi- 
nitions for s G [0, so}: 

V 2 j-i(s) = aj(s)aj(s) V 2j (s) = (3j(s)aj(s)* V (s) = j(s)a (s) 

where o!j(s) = a|(A(s)) for a smooth path of vectors A(s) such that < Xj(s) < 1 for s < sq. The 
50(2 J + 1) model is recovered at s = sq, where Aj = 1 and the parameters aj, Pj and 7 are given 
above. The commutation relations along the path can be read directly from those of the twisted 
CAR algebra above. Only in (4.9) do the parameters atj,/3j and 7 appear, namely 

(4.13) V 2j {s)V 2j - X {s) + \ 3 {sfV 23 -i{s)V 23 {s) = ^MfMvo(s) 2 . 

7 [S) 
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Using the commutation relations, the isometry condition for any s reads 



E (w 



\2 I |2 

Xj \aj\ 



which implies 
(4.14) 



E 



+ 7 



E 



and |/3j| 



Af | aj 



1, or equivalently 



IAf = l-Aj(l 



Of 



From (4.14) it is natural to interpret these coefficients as angles on the sphere S J , 

1 



with 



y(s) = — cos(s) and ctj(s) 



Vj 



sin s 



for all j, and cos(so) = (2J + 1) 1 / 2 . The PVBS algebra is recovered as s — > 0, since ctj(s) — > 0, 
while /3j(s) -»• 1 - A^O) 2 > 0. Since the left-right symmetry has not been broken along the path, 
the PVBS particles come in pairs with parameters Aj(0) for the left edge states and Aj(O) -1 for 
its right edge counterpart, see (4.11) and its adjoint relation. The singularity of the representation 
observed in the J = 1 case arises here too, as V2j(0) = 0. 

The intersection property for the spaces generated by the matrices {V a (s)}^L follows again from 
Lemma 2.3 with 

#(2) = {(00)} U {(a/3) : < a < /3 < 2J} . 

Hence, these MPS form the ground state spaces of Hamiltonians with nearest-neighbor interaction 
defined by the projection onto their orthogonal complement. 

In the next lemma we summarize the essential spectral properties of the transfer operator Ej(s), 
acting on -M 2 J "■ We note A > 0, resp. A > 0, whenever A%j > 0, resp. Aij > 0, for all i,j, and 
call these matrices nonnegative and positive, respectively. 

Lemma 4.2. We consider equivalently Ej(s) or its transpose for s > 0. 1 is a simple eigenvalue 
with a positive eigenvector. Moreover, all other eigenvalues have magnitude smaller than 1. 

Proof. We drop the s-dependence. Once again, it is convenient to consider the Ej for an arbitrary 
subset S C {1, . . . , J}, which is constructed from matrices {V^a-ii V^Jags and Vq , where is 
obtained from V a by keeping only the factors corresponding to indices in S. Let m = \S\. 

We decompose a general 2 m x 2 m -dimensional matrix A = A + (g> a + + A_ (g>a~ +Ap®P + Aq eg) Q 
to obtain 



E S j(A) 
(4.15) 



-X j Ef(A+) ® a + - XjEf(A^) ® ct" 



Ef (A P ) + \a j f(V ( ryA Q Vt 



sj ■ 



rS3 



P + 



WiyfyApvf + ix, 



Ef(A Q ) ®Q. 



First, the eigenvalue of Ej corresponding to an 'off-diagonal' eigenvector must be smaller than any 
eigenvalue of Ej j by a factor —Xj. But Ej J is also a principal submatrix of the diagonal part, which 
is nonnegative. Hence, its spectral radius must itself be smaller and we can restrict our attention 
to Ef tdia gonal- Its matrix representation has the following block form, 



M- 



m 2 D 



D 



\XA' Z M 



SJ 
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where D S3 is the diagonal matrix containing all possible products Y[i ± i k ^j \\ ' ' ' ^ for < A; < 
m — 1. If a, 6 index the blocks, it is immediate from 

K MS r] a , b = E < cl ---<- 1I 6 

Cl,...,C m _l=l 

that for any a, b, the submatrix M 5 ^ appears at least m — 1 times in at least one term of the sum. 
Since the other factors are diagonal with positive matrix elements, this implies that (M s ) m > if 
(M ) m ~ 1 > 0. Since the matrix depleted to just one particle reads 

we conclude by recursion that (M s ) m > 0. But this implies that the nonnegative matrix is 
irreducible, so that by Perron-Frobenius' theorem, its spectral radius is a simple eigenvalue with a 
positive eigenvector, and all other eigenvalues have a smaller magnitude. 

By construction, the identity matrix is an eigenvector of the full Ej for the eigenvalue 1, which 
implies that 1 is the spectral radius of Ej. This concludes the proof of the lemma for Ej. Its 
transpose shares its spectrum, and by nonnegativity and irreducibility, it must also have a positive 
eigenvector. □ 

Clearly, the proof collapses at s = as oy(0) = 0. In fact, the matrix M s becomes lower trian- 
gular and therefore reducible, and the eigenvectors fail to be positive. We recover the singularity 
of the path of the previous section: The positive diagonal left eigenvector pj of Ej, normalized to 
represent a density matrix, corresponds to a faithful state for s > 0, but not anymore at s = 0. 
Once again, this is necessary as the minimal representation of the PVBS model is one dimensional. 

The commutation relations allow for a simple description of the ground state space along the 
path of Hamiltonians. For each pair of particle indexed by j, the matrices V-ij and Vij-\ must 
alternate, separated by an arbitrary number of any other V a . If the chain is long enough, this 
produces A J different possibilities, depending on whether the first and the last matrix of each of 
the J types has en even or an odd index. Since the initial space of the matrix product map is -M 2 J ' i 
we have a complete d escr iption of the ground state spaces along the path. 



Part (i) of Lemma 3.8 holds for arbitrary J, as it requires only the spectral picture of Ej, the 
fact that 11 is the eigenvector for the largest eigenvalue 1, and the faithfulness of pj. In other words, 
the path is gapped for all s > 0. The proof of part (ii) would rely on the complete diagonalization 
of the transfer operator, which will not be carried out here. However, the key element is the 



'cost' of adding a pair above the PVBS ground state. From (4.13), it is precisely ajfij/'jj which 

is proportional to tan(s)/(\/«/ cos(s)) for all j, indicating that the perturbative argument can be 
carried out in the general case. This would imply that the gap does not close as s — > 0, and the 



generalization of Theorem 3.1 For any J £ N, the SO(2J + l)-invariant model belongs to the 
PVBS phase with = tir = J. 

The paths constructed here and the necessary singularity encountered at the PVBS point raise the 
more general question of the continuity of the spectral gap of an interaction exposing a continuous 
family of matrix product states, in the spirit of [21J. The fundamental difficulty lies in the non 
invertibility of the eigenvectors of E at the PVBS point, and therefore in the collapse of standard 
estimates involving p _1 and e . Such 'stability' results are notoriously rare (see |22| for a concrete 
example) but we do believe that continuity holds more generally than in the simple case presented 
here. This is subject to current investigation. 
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